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Knowledge of the derivative of the topological susceptibility at zero momentum is important for 
assessing the validity of the Witten-Veneziano formula for the 77' mass, and likewise for the res- 
olution of the EMC proton spin problem. We investigate the momentum dependence of the topo- 
logical susceptibility and its derivative at zero momentum using overlap fermions in quenched 
lattice QCD simulations. We expose the role of the low-lying Dirac eigenmodes for the topologi- 
cal charge density, and find a negative value for the derivative. While the sign of the derivative is 
consistent with the QCD sum rule for pure Yang-Mills theory, the absolute value is overestimated 
if the contribution from higher eigenmodes is ignored. 
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1. Introduction 

The derivative at zero momentum of the topological susceptibility x{^^) iri momentum space, 
;^'(0), is important to control the validity of the Witten-Veneziano formula for the r\' mass [|l], 
The derivation of this formula presupposes that the momentum dependence of the topological 
susceptibility x{^^) is moderate from = 0to m^,, i.e. \x'm^i\ ^ X- Moreover, the large-A'^c hmit 
is implied. The derivative is also helpful to analyze the EMC proton spin problem [||] (see [Q] for a 
review). Therefore it is of great interest to investigate x'i'^) starting from QCD. 

Based on various kinds of QCD sum rules, the following estimates of ;^'(0) have been reported. 
For pure gauge (Yang Mills) theory, x'{'^)\ym = — (7±3 MeV)^ and with account of dynamical 
quarks, 7^(0) « 23.2 ± 2.4 MeV [§, x'{0) = (2.3 ± 0.6) x 10"^ GeV^ [^, ^/x^ = 26.4 ± 
4.1 MeV ix' ^ 0.7 x 10"^^ GeV^) [|], and x'{0) « 1.82 x lO^^ GeV^ 

On the other hand, there are few lattice QCD simulations, which resulted in the following 
values. In quenched lattice theory, x'(0)|su(2) = -(9.84±0.91 MeV)^ [Q, ;t:'(0)|su(3) = -(13 ± 
16 MeV)2 [|ll|], and full lattice QCD with staggered fermions, a/zW = 19 ±4 MeV [^]. In these 
simulations, the topological charge density has been defined using the (unimproved) field strength 
tensor constructed in terms of plaquette variables, and the cooling method has been applied to 
eliminate the ultraviolet (UV) noise. Here, one may notice that the sign of the estimated by 
the QCD sum rules is negative in pure Yang-Mills theory while it is positive in full QCD, and the 
lattice results seem to support this tendency within numerical errors. 

In this paper, we are going to use a different approach in order to further investigate the deriva- 
tive of the topological susceptibility at zero momentum using lattice QCD simulations. We shall 
use overlap fermions [|^, |l^] to define the topological charge density [jl5|]. This approach is ex- 
tremely useful to clarify the topological structure of the QCD vacuum, since it preserves exact 
chiral symmetry on the lattice JI^ and satisfies the index theorem [[l7|]. In this approach, it is also 
possible to expose the role of the low-lying Dirac eigenmodes for the topological structure of the 
QCD vacuum [|l|]. 



2. Topological charge density and susceptibility from overlap fermions 



In order to define the topological charge density q{x) on the lattice with the lattice spacing a, 



we employ the massless overlap Dirac operator [13, 14] defined by 



D=P(l + 



X 



X 



a 



(2.1) 



where Dyv is the Wilson-Dirac operator. We set p = 1 .4, a value identified for the lattices in use as 



an optimal choice. Details of our implementation are described in ||18|]. Overlap fermions possess 
exact chiral symmetry on the lattice JT^ ] and provide ?i_ +«+ exact zero modes, = 0, with 
(«+) being the number of modes with negative (positive) chirality: 75 V'^ = — Wn ^^i^ YiWn — ~^¥n- 
The index is given by Q = — [17]. The non-zero modes with eigenvalue A, D^x = ^Y^^ 
occur in complex conjugate pairs X and X * and satisfy Y^TcP^six) = Y.^^ Tr Yx i^) 75 ¥x (■^) = 0, where 
Tr should be regarded as the sum over color and spinor indices. 
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Then, the topological charge density q{x) is given by [ p^ ] 

^(^) = _Tr[75(l-^D(x,^)) 



(2.2) 



which satisfies the index theorem with Q = Y.x^{^) ^ ^- The contribution of the low-lying Dirac 



eigenmodes to q{x) can be exposed by applying the eigenmode expansion [19, 20, 21], 



'?(^;Acut) 



1-^ IpasW 



(2.3) 



where the cut-off Acut implies a kind of UV filtering, and pxsi^) = Wx{^)7iWx{^) the local 
chirality of the mode corresponding to the eigenvalue A . Note that the UV filtering maintains the 
index theorem independently of the cut-off such that Q = Y^x^i^'^^caij- This is because the index 
is computed only from the zero modes. Actually, it is observed that all zero modes of a lattice 
configuration occur with the same chirality. Truncating the expansion at Acut acts like an UV filter 
by removing certain short distance fluctuations from the local density q{x). 

In the continuum limit, the momentum dependent topological susceptibility is defined as 



(2.4) 



where Cq{x) = {T {q{x)q{0))) is the two-point correlation function of the topological charge density 
^(x). Eq. (2^) can be expanded with respect to and one obtains the expression of the derivative 



of the topological susceptibility at zero momentum 



x'io) 



dx{k^) 



dk^ 



d^X X^Cq{x) 



(2.5) 



The lattice version of Eqs. (2.4) and (2.5) reads 



;t(P) = £.'%(x), 

xev 



(2.6) 



respectively, where 



° xev 



xeV 



(2.7) 



(2.8) 



As we consider the lattice volume V = L^T and impose periodic boundary conditions in all direc- 
tions, the momentum is discrete as 



2 . , dJCn 



with being integer; ,?i2>'^3) runs from to (L/a) — 1 and ?i4 from to {T /a) — 1. 
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Table 1: Details of the ensemble used in this study. The Liischer-Weisz gauge action is used, j3 denotes the 
inverse coupling squared, a the lattice spacing (determined from the pion decay constant), and V =L?T the 
lattice volume. ^(0) is the topological susceptibility. 





a [fm] 


{L/af{T/a) 


(fm^) 


# of config. 


[GeV^] 


8.45 


0.105 


16^32 


(15.9) 


267 


8.2(7) X 10^4 


8.45 


0.105 


12^24 


(5.0) 


116 


9.1(13) X 10-4 


8.10 


0.142 


12^24 


(16.9) 


254 


8.6(8) X 10-4 


8.00 


0.157 


16^32 


(79.6) 


444 


9.1(6) xlO-4 



3. Numerical results 

We used several ensembles of zero-temperature quenched configurations generated by means 
of the Liischer-Weisz gauge action [22] (see Table |I|). This action is suitable for topological studies 



since dislocations are greatly suppressed. 

We plot xiP') as measured at j3 = 8.45 on the 16^32 lattice in Fig [l] (left), and as measured 
at j8 = 8.10 on the 12^^24 lattice in Fig. [I] (right), for various cut-off values of the (imaginary part 
of the) eigenvalue, Acut = 0.2 — 0.8 GeV. In these analyses, the eigenvalue A of Eq.( |2.1[ ) is 
always replaced by Ajmp = (1 — aA/2p)-'A. These are the eigenvalues of the improved massless 
overlap operator Dimp [p^]. The improvement projects the eigenvalues of D stereographically onto 
the imaginary axis. We find that the largest value of ;^(P) is at P = for each cut-off value of 
the eigenvalues Acut = 0.2 — 0.8 GeV and that the functions ;t(P) are monotonously decreasing 
as functions of the momentum for all cut-off values. The qualitative behavior of xi}?') does not 
depend on the lattice spacing. 



The result of the x'{^) according to Eq. (2/7) is plotted in Fig. g. We do not see a significant 



volume dependence at j8 = 8.45 for which two lattice volumes are available. Since the physical 
volume of the 12^24 lattice at j8 = 8.45 is the smallest among our lattice ensembles (see Table [I]), 
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Figure 1: The momentum dependence of the topological susceptibility at )3 — 8.45 on the 16^32 lattice 
(left). There are two data points at some momenta, which are due to the violation of rotational symmetry, 
and are expected to converge to one point in the continuum limit. The same as the left figure, but at j3 =8.10 
on the 12^24 lattice (right) . 
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0.002 



we expect that the finite volume effect will not spoil our observation. The physical scales of the 

cut-off values at different j8 values are chosen approximately corresponding to each other, so as to 

examine the lattice spacing dependence of x'(0)- In the result we find a reasonable scaling behavior 

of x'{^) corresponding cut-off values Acut, apart from the data of the coarsest lattice at j8 = 8.00. 

In any case, we find an interesting behavior: as the cut-off increases, |;t'(0)| becomes smaller. 

In other words, |;t'(0)| takes on the largest values if only the lowest eigenmodes are included in the 

fluctuations of the topological density. If one, considering the Witten-Veneziano formula, ignores 

the necessary condition of \x'{Q)m'^,\ ^ ;t(0), it seems that the topological susceptibility at zero 

momentum ;^(0) is the only important parameter for the T]' mass. Since ;^(0) depends only on the 

(number and chirality of the ) zero modes, one might erroneously conclude that the zero modes are 

sufficient to describe the topological structure of the QCD vacuum. However, the behavior of ;t'(0) 

indicates that fluctuation of the topological charge density at all scales are necessary to warrant that 

\x'{^)^\'\ < Z(0) can be satisfied. 

For instance, taking a typical value at 

j8 = 8.10, we have ;t;'(0) ^ -0.001 GeV^ 
Ct;'(0) ^ -(32 MeV)2) with ;t;(0) ^ 8.6 x 
lO^-* GeV^, we estimate z(0)/|z'(0)l ~ 
0.86 GeV2 = (930 MdWf, which is com- 
parable with m^i = 958 MeV. This in 
turn means that not only the lowest eigen- 
modes (not to mention only the zero 
modes) but also the higher eigenmodes 
have to play an important role for the 
topological structure of the QCD vacuum. 

It is then quite intriguing to inves- 
tigate which value the derivative 
approaches if more and more contribu- 
tions from higher eigenmodes (possibly 
all eigenmodes) are included. We in- 
vestigate this for the case of j8 = 8.10 
on the 12^24 lattice based on the defini- 
tion of the full topological charge density 
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Figure 2: The slope ;i;'(0) as a function of the cut-off 

Acut- 



5x10" 



in Eq. (2.2), which was evaluated with- 



out mode expansion. We find the value 
of flie derivative ;t'(0) = -2.4(16) x 
lO^^* GeV^. However, we note that the 
number of configurations used in this 
analysis is only 53 [|l^] and that the 
topological susceptibility is overestimated 
for this subsample as ;^(0) = 1.3(3) x 
10^^ GeV^. The limited number of The 
small number of configurations is dictated 
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Figure 3: The same as in Fig. 1, but for the ensemble 
with j3 = 8.10 on the 12^24 lattice, when the zero modes 
and the lowest eigenmodes are removed from the eval- 
uation of Xv^^)- Only eigenmodes between 0.4 < A < 
0.8 GeV are taken into account. 



by the high cost of evaluating Eq. (2.2) us 



5 



Momentum dependence of the topological susceptibility with overlap fermions 



Yoshiaki Koma 



ing the overlap Dirac operator. Therefore, the value of ;^'(0) based on the all-scale topological 
density (formally Acut — ?• °° ) cannot be considered as conclusive. 

Next we examine the role of the zero and lowest modes for xi^^)- I^i Fig. |3|, we plot the 
X(k^) for the ensemble with jS = 8.10 on the 12^24 lattice, where only the eigenmodes between 
0.4 < A < 0.8 GeV are used to define the topological charge density, i.e., the zero modes and very 
low -lying modes are removed. We find that the slope at P = turns to be positive. ;t;(P) increases 
from zero at ^ = (due to the subtraction of zero modes from the topological density) to a certain 
maximum value (near to P 1 GeV^ ) and then decreases. 

This leads us to the speculation that - if the zero modes are dynamically suppressed and the 
relative weight of sectors with smaller \Q\ will increase in the average - the sign of the derivative 
X'iO) may change from negative to positive. This may be the case in full QCD in the chiral limit. 

4. Summary 

Using the overlap fermion formalism, we have investigated the momentum dependence of the 
topological susceptibility and its slope at zero momentum, ;t'(0). Overlap fermions preserve exact 
chiral symmetry on the lattice and possess exact zero modes, which allow us to unambiguously 
compute the index Q of vacuum configurations. We have found that x'i'^) depends on the number 
of eigenmodes determining the resolution in the definition of the topological charge density. The 
more the cut-off applied to the eigenvalues of the Dirac operator increases, the more the absolute 
value of the slope |x'(0)| decreases. In other words, |x'(0)| with only the lowest-lying eigenmodes 
included turns out too large to guarantee that the technical assumption underlying the Witten- 
Veneziano formula is fulfilled. From this point of view, fluctuation of the topological charge density 
at all scales, represented by higher eigenmodes, should also be taken into account. Therefore it is 
not correct to argue that only zero modes and low-lying eigenmodes are relevant for the topological 
structure of the QCD vacuum. 

As demonstrated in |[T^], the topological charge density possesses global sign coherent struc- 
tures, which get increasingly tangled as more and more eigenmodes are included, and the all-scale 
topological charge density has a lower-dimensional, laminar structure, together with a lumpy struc- 
ture inside the sign coherent regions. After all, it is found to possess a multifractal structure. This 
complicated structure of the QCD vacuum is responsible for the behavior of x'i^)- 

The numerical calculations have been performed at NIC (Jiilich) and HLRN (Berlin) as well as 
at DESY (Zeuthen). We thank all institutes for support. This work was supported in part by Japan 
Society for the Promotion of Science (JSPS) and German Research Foundation (DFG), Japan- 
German Joint Research Project 2008-2009. Y.K. was partially supported by the Ministry of Edu- 
cation, Science, Sports and Culture, Japan, Grant-in-Aid for Encouragement of Young Scientists 
(B), No.20740149. M.K. was supported by Japan Society for the Promotion of Science (JSPS), 
Grant-in-Aid for JSPS Fellows (20-40152). 
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